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v  Abstract 

The  steady-state  distribution  function  is  obtained  for  electrons 
initially  emitted  from  a  point  source  into  a  neutral  gas  and  which  subse¬ 
quently  drift  under  the  influence  of  a  uniform  dc  electric  field  while 
undergoing  elastic  collisions  with  the  gas  atoms.  The  usual  approxi¬ 
mations  are  retained  of  regarding  the  distribution  function  as  almost 
spherical  in  velocity  space  and  of  regarding  the  fractional  energy  gain 
or  loss  by  an  electron  upon  collision  as  small.  However,  the  terms  in 
the  Boltzmann  transport  equation  involving  spatial  derivatives  of  the 
distribution,  which  are  usually  assumed  small  in  comparison  to  the 
field  and  collision  terms,  are  treated  exactly.  The  distribution  function 
is  given  as  a  sum  of  energy  modes,  each  of  which  decay  with  distance 
from  the  source.  The  lowest  of  these  modes  is  the  far-distant  distribution, 

't 

while  the  higher  ones,  which  decrease  more  rapidly  with  distance,  describe 
the  decay  of  the  initial  source  energy  distribution.  The  complete  distri¬ 
bution  is  obtained  in  terms  of  known  functions  in  the  case  of  an  energy 
independent  collision  frequency,,  whereas  in  the  energy  independent 
cross  section  case,  only  the  lowest  mode  is  obtained.  The  far-distant 
part  of  the  distribution  function  is  compared  with  the  usual  approximate 
expression  which  is  obtained  when  the  gradient  terms  are  considered  small 
and  which  is  expressed  as  the  density  times  a  normalized  energy  function. 
It  is  shown,  that  when  the  gradient  terms  are  correctly  considered,  the 
far-distant  distribution  in  energy  becomes  position  dependent.  Furthermore, 
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Abstract  (Cont. ) 

the  deviation  from  the  approximate  theory  becomes  larger,  the  further 
the  electrons  are  off  the  geometrical  axis.  This  position  dependence 
is  most  important  when  the  electron  energy  is  large  in  comparison  to 
thermal  energies.  The  interpretation  of  Townsend  method  for  the 
determination,  of  the  ratio  of  the  diffusion  coefficient  to  the  mobility, 
D/jjl,  is  re-examined  on  the  basis  of  this  more  exact  theory.  It  is 
shown  that  the  error  in  D  j |a  that  results  from  using  the  conventional 
interpretation  of  this  method  under  typical  experimental  conditions  is 
never  more  than  about  20%. 
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I.  INTRODUCTION 

It  is  usually  assumed  that  electrons  which  drift  and  diffuse 

through  a  gas  under  the  influence  of  both  uniform  electric  fields  and 

electron  density  gradients  have  a  distribution  in  energy  that  is 

independent  of  position.  This  means  that  the  distribution  is  assumed 

to  be  unaffected  by  the  presence  of  gradients  in  the  electron  density  and 

is  taken  to  depend  only  on  the  field  strength  and,  of  course,  on  the 

pressure  and  variety  of  the  gas.  A  direct  consequence  of  this  is  that 

the  electrons  can  be  characterized  by  a  diffusion  coefficient,  D,  and 

a  mobility,  p,  which  are  independent  of  position.  It  is  this  aspect 

that  has  been  the  basis  for  the  interpretation  of  many  experiments 

concerned  with  the  transport  properties  of  electrons  in  gases,  such  as 

/  1,  2 

the  Townsend  type  experiment  for  the  measurement  of  D/  p  ’  and  the 

1  3 

time  of  flight  measurements  of  p.  ’  In  turn,  the  analysis  of  the 

measured  transport  coefficients  from  these  experiments  in  terms  of 

electron-atom  collision  cross  sections  have  also  been  dependent  on  this 
-1,4 

assumption. 

The  theoretical  justification  of  this  assumption  must  come 
from  the  solution  of  the  Boltzmann  transport  equation  that  is  appropriate 
to  electrons  under  the  influence  of  both  electric  fields  and  electron 
density  gradients.  Allis  and  Allen5  have  derived  the  basic  equations  for 
electrons  under  these  conditions.  These  authors  did  indicate  certain 
formal  aspects  of  the  general  solution  as  well  as  pointing  out, 
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qualitatively,  the  approximate  nature  of  the  conventional  assumption; 
however,  they  did  not  discuss  any  specific  case  fully. 

It  is  the  purpose  of  the  present  study  to  (1)  investigate  the 
conditions  under  which  the  effect  of  electron  density  gradients  on  the 
distribution  in  energy  can  be  neglected  and  (2)  to  obtain  solutions  of  the 
Boltzmann  equation  which  will  demonstrate  the  specific  effects  intro¬ 
duced  by  the  gradients.  The  geometry  used  in  this  study,  which  is 
one  of  the  simplest  that  can  serve  to  illustrate  these  effects,  is  the 
point  source  of  electrons  in  an  infinite  uniform  field  region.  Also 
for  reasons  of  simplicity,  the  collisions  between  the  electrons  and 
gas  atoms  are  taken  to  be  elastic.  In  Section  II  the  Boltzmann  equation 
for  the  case  of  an  energy  independent  collision  frequency  is  given  and 
the  approximations  involved  in  assuming  the  distribution  in  energy 
to  be  position  independent  are  examined.  Also  a  criterion  is  developed 
for  the  conditions  under  which  the  effect  of  electron  gradients  can  be 
neglected.  In  Section  III  the  Boltzmann  equation,  as  given  in  Section  II, 
is  solved.  The  resulting  solution  yields  not  only  the  limiting  distri¬ 
bution  at  far  distances  from  the  source  but  also  the  higher  modes  which 
describe  the  decay  of  the  initial  distribution  into  this  far-distant  part. 
The  far-distant  distribution  is  then  compared  with  the  usual  position 
independent  distribution.  In  Section  IV  the  Boltzmann  equation  for  the 
case  of  an  energy  independent  cross  section  is  given  and  an  approxi¬ 
mate  expression  is  obtained  for  the  distribution  at  far  distance  from  the 
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source.  This  expression  is  compared  with  that  obtained  for  constant 
collision  frequency.  In  Section  V  the  Townsend  D/ p.  experiment  is 
examined  on  the  basis  of  this  more  exact  theory  and  the  errors  that 
result:  from  using  the  conventional  interpretation  of  this  experiment 
are  discussed. 


II.  GENERAL  CONSIDERATIONS 

The  average  properties  of  electrons  moving  through  a  gas, 
e.  g.  ,  electron  density,  current  density,  mean  energy,  etc.  ,  can  be 
predicted  once  the  electron  distribution  function,  f(r,  v),  is  known. 

<w  /w 

The  significance  of  this  function  is  that  f(r,  v)  dr  dv  denotes  the 

/V Al  |VM  /w 

number  of  electrons  at  position  r;  in  dr  and  with  velocity  v  in  the  range 
dy.  The  distribution  function  in  turn  satisfies  an  equation  of  continuity 
in  position  and  velocity  space,  i.  e. ,  the  Boltzmann  transport  equation. 
This  equation  describes  the  balance  that  must  exist  in  steady  state 
between  the  rate  at  which  electrons  enter  and  leave  a  given  element 
of  volume,  dr;  dv;  in  velocity  and  position  space.  The  flow  in  position 
space  results  from  the  velocity  of  the  electrons  while  in  velocity 
space  it  results  from  their  acceleration  due  both  to  collisions  with 
the  gas  atoms  and  to  the  applied  field. 

There  are  several  approximations  commonly  made  in  order 
to  simplify  the  integral-differential  Boltzmann  equation  when  applied 
to  electrons.  The  first  is  that  the  distribution  function  is  almost 
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spherically  symmetric  in  velocity  space  and  therefore  can  be  adequately 
represented  by  the  first  two  terms  of  an  expansion  in  spherical  har¬ 
monics  involving  the  direction  of  the  velocity.  That  is,  f(r,  v)  can  be 

,  6 
written  as 


f{r,  v)  ='  f°(r,  v)  +  fV.  v)  •  (v) 

*/✓«/■  Ml  AV 


The  second  is  that  the  fractional  energy  gain  or  loss  by  an  electron 
upon  colliding  with  a  gas  atom  is  small.  This  is  justified  in  the  case 
of  elastic  collisions,  to  which  the  present  paper  is  restricted,  because 
of  the  small  electron  to  atom  mass  ratio.  With  these  approximations 
the  Boltzmann  equation  reduces  to  two  partial  differential  equations. 

The  present  discussion  will,  in  addition,  be  restricted  to  uniform  dc 
electric  fields  and  to  a  constant  collision  frequency  gas.  While  this 
latter  restriction  will  be  relaxed  further  on  to  include  the  case  of  con¬ 
stant  cross  section,  for  the  present  the  case  of  constant  collision  fre¬ 
quency  can  serve  best  to  illustrate  the  important  features  of  the  problem. 

When  the  above  approximations  and  restrictions  are  taken 

o  16 

into  account  the  equations  that  result  for  f  and  f  are 


2ro.v  J3 _  rV 2  fr°  x  v  t 
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Here  the  electric  field  ^is  given  by  E^=  -kE,  where  k  is  the  unit 

vector  in  the  z  direction,  v  is  the  momentum  transfer  collision  frequency, 

£  denotes  the  kinetic  energy  of  the  electrons,  m  and  M  are  the  mass 

of  the  electron  and  the  atom  respectively,  T  is  the  gas  temperature, 

and  S(r,  €)  is  the  electron  source  term, 

v. 

Qualitatively  these  equations  can  be  explained  as  follows: 

We  have  assumed,  that  f(r,  \j)  can  be  represented  by  the  term,  f°,  that 
is  spherically  symmetric  in  velocity  plus  the  small  non-spherical 
term,  •  v.  Therefore  the  Boltzmann  equation,  which  balances 
the  rates  at  which  electrons  enter  and  leave  dr  dv,  also  breaks  into 
two  parts,  the  first,  Eq.  (1),  which  balances  the  spherical  rates  and 
the  second,  Eq.  (2),  which  balances  the  non-spherical  rates.  The 
first  term  in  (2),  which  represents  the  effect  of  collisions  in  reducing 
the  asymmetry  in  f(rv),  is  balanced  by  the  second  and  third  terms 
which  represent,  respectively,  the  effect  of  diffusion  and  drift  in 
increasing  the  asymmetry.  The  first  two  terms  in  (1)  reflect,  respec¬ 
tively,  the  fact  that  electrons  can  lose  and  gain  energy  as  a  result  of 
collisions  with  the  gas  atoms.  The  third  term  represents  the  net 
flow  of  electrons  into  d^with  energy  £  that  occurs  when  the  asym¬ 
metrical  part  of  f  varies  with  position.  The  fourth  term  reflects  the 
fact  that  the  electrons  can  gain  energy  from  the  field  and  that  this 
occurs  only  through  the  asymmetrical  part  of  f. 
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The  equation  that  f  must  satisfy  is  obtained  by  substituting 


f1 


f^  from  (2)  into  (1)  and  is 
2mv  3 


M£' 


■1/2  Be 


Bf° 

B£ 
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MeE 


{-) 

\mv  / 


2Jt£ 

3  2 


I 
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,  2eE£  S2f°  _2jE  2  o 
3mv„  3mv  v* 


=  -S(r,6)  . 


(3) 


O 


When  f°  is  independent  of  position,  i.  e.  ,  when  the  electron  density  is 
uniform,  the  solution  of  (3)  with  S  =  0  is  given  by 


Const  x  exp 


(4) 


Under  these  conditions  the  balance,  as  represented  by  the  first  three 
terms  of  (3),  is  between  the  electrons  losing  energy  from  collisions  and 
gaining  energy  from  collision  and  from  the  field.  When  f°  depends  on 
position,  the  additional  terms  in  (3)  that  involve  the  spatial  gradient  of 
f°  appear.  If  these  terms  are  assumed  to  be  small  in  comparison  to  the 
collision  terms  and  to  the  field  term,  then  an  approximate  solution  to  (3) 
can  be  expressed  as  f°  ~  n  (|j)  F(£  ),  where  n  (£)  is  the  electron  density, 
whose  functional  form  is  as  yet  undetermined,  and  F(6  ),  which  repre¬ 
sents  the  distribution  in  energy,  is  given  by  (4).  The  equation  that  n(r) 

1  /  2 

must  satisfy  is  obtained  by  multiplying  (3)  by  £  d£  and  integrating 
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over  the  complete  energy  range.  The  resulting  equation,  which  is  the 
familiar  continuity  of  current  equation,  can  be  written,  with  S  =  0,  as 


(D^n  -  tiEkn)  =  <3, 


(5) 


where  D  and  |x  are  given  by 


“•3  e>*  J 


.3/4 


F(€)  d6 


and 


81Te  f  £3/2  dF(€  ) 

3m2  J  v  d£ 

o 


d£  . 


A  simple  criterion  can  be  obtained  for  conditions  under  which 
the  gradient  terms  can  be  neglected  in  comparison  to  the  field  term  or  the 
collision  terms.  If  the  approximate  form  for  f°  is  assumed,  i.  e. , 

f°  -  n(r^  exp  € /(kT  +  1/4)  ^  ) 

where  B  =  (3/M)  (mv/eE)^  arid  this  is  substituted  into  (3)  along  with  the 
2 

equation  for  n  given  by  (5),  then  the  relative  magnitude  of  the  various 
terms  can  be  compared.  For  this  case  of  constant  collision  frequency  the 
last  two  terms  on  the  left  hand  side  of  (3)  cancel.  The  ratio,  of  the 
remaining  gradient  term  to  either  the  energy  loss  collision  term  or  to 
the  sum  of  the  collision  and  field  term  representing  energy  gain  can  be 


expressed  as 
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•  (7> 

In  this  equation  €av  denotes  the  average  electron  energy,  (3/2)(kT+l/B), 
is  the  thermal  energy  of  the  gas,  (3/2)  kT,  and  D/ p.  =  (kT+l/B)/e  . 

We  see  that  if  € >eT  (high  field  limit),  the  gradient  terms  can  be 
neglected  when  the  diffusion  current  in  the  field  direction,  DOn/fez), 
is  small  in  comparison  to  the  drift  current,  uEn.  However  as  C 

av 

approaches  £  (low  field  limit)  the  gradient  terms  become  less  and  less 
important  for  a  given  ratio  of  diffusion  to  drift  current. 

The  case  of  a  point  source  in  an  infinite  uniform  field  region 
can  serve  to  illustrate  these  points  for  a  specific  geometry.  The  density 
for  a  point  source  is  given  by  ^ 


2  2  -V2 

n(r)cC  (z  +  p  ) 


where  z  is  the  distance  from  the  point  source  along  the  field  direction 
and  p  is  the  cylindrical  radius.  When  the  density  gradient  to  density 
ratio  is  obtained  from  this  expression  and  is  substituted  into  (7),  the 
ratio  of  terms  in  the  high  field  limit  (i/b>  kT)  becomes 


cos  8  - 


4  C’ 

w  av 

3eEz 


2 

cos 


where  9  is  the  polar  angle  from  the  source.  When  the  factor  4  £ 
is  small  and  the  point  of  interest  is  near  the  axis  (cos  9*1),  the  gradient 


.av/ 


/3eEz 
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term  can  be  neglected.  However,  for  positions  sufficiently  off  the 

axis  (cos  0  <  1)  and  again  with  4  C  / 3eEz  ^1,  the  gradient  term 

a  v 

cannot  be  neglected  under  any  circumstance. 

The  extension  of  this  type  of  argument  to  a  constant  cross 
section  gas  is  straightforward  and  probably  the  above  conclusions  are 
a  reasonable  guide  for  gases  with  a  more  complicated  energy  dependent 
collision  frequency. 


III.  CONSTANT  COLLISION  FREQUENCY 


A.  General  Solution 


The  equation  for  f  ,  which  is  given  by  (3),  can  be  rearranged 


to  read 


2,o 


_i_  s>zf' 

eEB  : 


(eE)2B 


T  f' 


2fo  =  _  Mg(fcg  l 
2mv  £ 


where  B,  as  defined  previously,  is  ’  (^/mJ^j/^E)2.  It  is  clear  that 
because  of  the  mixed  derivatives  in  £  and  z  the  equation  is  not  separable 
in  these  variables.  However,  by  changing  to  new  independent  variables, 
that  are,  for  convenience,,  made  dimensionless,  the  equation  can  be  put 

g 

into  a  separable  form.  The  new  variables  are  defined  as  follows, 
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x  = 


B€ 


1  +C* 


* 


=  x  -  eEBz  , 


R  =  eEBp  , 


and  with  Ok  =  kTB.  Then  the  differential  equation  becomes 


Sf° 


2f°  ^ 

-2  +  (1+o()- 


9r 


MS 

2mvx 


(9) 


where  has  been  expressed  in  cylindrical  coordinates  since  f°  is  to 
be  calculated  under  conditions  of  cylindrical  symmetry.  It  should  be 
pointed  out  that  in  the  high  field  limit,  i.  e.  ,  when  o«y.  the  new 
variable  is  simply  proportional  to  the  total  electron  energy,  6  -eEz. 
Since  the  equation  separates  in  the  variables  x  and  ^  and  not  in  6  and 
z,  it  would  be  a  formidable  task  to  satisfy  boundary  condition  on  a  z-plane. 
However,  what  can  be  obtained  in  a  straightforward  fashion  are  solutions 
for  the  case  of  volume  sources  in  an  unbounded  region.  Another  case 
that  can  be  simply  treated,  but  will  not  be  detailed  in  the  present  paper, 
is  for  volume  sources  in  the  presence  of  cylindrical  boundaries  parallel 
to  the  z  axis  and  on  which  f°  vanishes,  e.  g. ,  a  point  source  on  the  axis 
of  a  right  cylinder  with  f°  vanishing  along  the  inside  surface.  However, 
the  simplest  geometry  that  can  be  used  to  illustrate  the  features  of  the 
distribution  when  the  gradient  terms  are  correctly  considered  is  the 


point  source  emitting  monoenergetic  electrons  in  an  unbounded  region. 
Therefore,  Eq.  (9)  is  to  be  solved  with  S  given,  in  terms  of  £  ,  z  and 
p,  as 


S(€,  z,  p)  =  S(z) 


S  (p) 


6(6  -  6) 


^P  (4Tf/m)(2  6/m)1//2 


(10) 


This  term  represents,  one  electron  per  second  being  emitted  with  energy 
from  the  point  z  =  p  =  o. 

The  homogeneous  form  of  Eq.  (9)  can  be  separated  into  the 
three  ordinary  differential  equations, 


x_3/2  ^  [x3/2  (F(x)  +i5dfSi)| +  P  F W  =  o  • 

Ri  (R^)  +  K2p(R)  '  0  ’ 

+  ^  + o  -  o  • 


(11) 

(12) 

(13) 


The  solutions  for  (12)  and  (13)  that  are  of  interest  can  be  written  down 
immediately  as 


P(R)  =  JQ  (KR), 


and  Z(Sj’)  =■  exp 


1  ±Cl+4Q(fK2  (1+OQ  +  b)]1^2 

2  04 


where  Jq  (KR)  is  the  zero  order  Bessel  function.  In  order  to  put  (11) 


in  a  recognizable  form  let 
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_JLul 


and 


P  = 


F  =  h(x)  exp 


L-<^  •]  ■ 


Then  (11)  becomes 


x 


d2h 


+  (3/2  ~vx)  “  -3/4  (fy  -  1)  h  =  0  . 


dx 


The  solution  of  interest  is  regular  at  the  origin  and  is 


9 

where  F(a(  c  (z)  is  the  confluent  hypergeometric  function.  Now  Eq.  (11) 
is  an  eigenvalue  equation  for  P  (or  y).  It  is  easily  shown  that  the 
corresponding  eigenfunctions  are  orthogonal  with  respect  to  the  density 
function  r(x)  =  x2/^2eX,  i.  e.  , 


x  '  2eX  F  F  ('  dx  =  o  if  y  f  V 

The  eigenvalues  and  eigenfunctions  are  found  by  selecting 

the  set  of  y 's  which  made  the  solution  of  (11)  quadratically  integrable 

and  orthogonal  with  respect  to  the  density  function.  By  inspection  of 

9 

the  asymptotic  behavior  of  the  confluent  hypergeometric  function  the 
allowed  spectrum  for  y  and  the  corresponding  eigenfunctions  can  be 
found.  A  part  of  the  set  is  discrete  with 
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5e 


M  =n7  exp[- 


X 


4/2V> 


(1/2)  9 

where  'is  a  Laguerre  polynomial, 


\  1  +  (4/3)j?  ’ 

Jl  =  0,  1,  2  .  .  .  . 


00  (positive  integers)  , 


9 

id  with  the  normalization  constant  given  by 

,T2  (1  +  (4/3 )Jl)5//2  (2 Jl  +  3/2)  T3  (3/2 

nji  j?  : 


The  rest  of  the  set  is  continuous  with 


CO 


F  (3/4 


where  Cwis  a  normalization  constant,  y  =iCoand  W  is  a  continuous 
variable  in  the  range  of  0  -►  +  00.  It  will  be  assumed  that  this  is  a 
complete  set. 

The  solutions  to  Eqs.  (11),  (12),  and  (13),  as  given  above,  will 
now  be  used  to  build  up  the  solution  to  Eq.  (9).  It  is  convenient  first, 
to.  find  the  Green's  function  for  (9),  i.  e.  ,  to  find  the  function 
G(x^R  j  x^R)  that  satisfies  the  equation 


x 


■3/2  9_ 

3x 


isL  (  r2*g.  \ 

r3r  \Br) 


-  £>(x  -  x)  S  )  S(R  -  R)  . 


(14) 
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and  then  to  use  the  expression 

f°  (x, ^ ,  R)  =  Jg(x, ^ ,  R | x ^R)  j^mv'f  * ^  dx  dR  (15) 

to  obtain  the  f°  that  satisfies  (9)  for  a  specific  source  term.  Now  the 
function  G  can  be  expanded  in  terms  of  the  energy  functions,  ,  as 

G  =«sEl  (^,  R,  R,x):F^(x), 

where  the  indicated  summation  is  to  be  taken  as  a  sum  over  the  discrete 
spectrum  of  y  plus  an  integral  over  the  continuous  part.  When  this 


expansion  is  substituted  into  (14)  and  the  result  is  multiplied  by  r(x)F^/(x) 
and  integrated  over  x,  the  equation  for  A^  is  found  to  be 


A 

-PA  +  + 

y  d%2 


»+«>  kk  (?w)m  ■  r®\  ® 


(16) 


Now  A  can  be  expressed  as 
co 


A  =  U  (^,  R^vx)  J  (KR)  KdK 


V  J  Y 

o 

When  this  is  substituted  into  (16)  and  the  result  multiplied  by  RJq(!!CK) 
and  integrated  over  the  complete  range  of  R,  the  equation  for  that 


results  is 


3u  c?u 

"+o<-2-X=-r 


-(p  +  (i+«)K)uv+-g? 


r®ys>  X<kR)  S(yt() 


The  solution  to  this  equation  is 
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exp 

=  R  r(x)  Mjffi — 


(V) 


exp- 


[l  +  4  (p  +  (l+0<)K 
Therefore,  G(x,^  ,  R  |x,  ^  ,  R)  is  given  by 


fi+4^([3  +  (l+OOK2D 
2<X 

W 


1/2 


G  =  R  exp 


[<!  -?>/20<]^ 

Yl  +  4  0(0  +  (1  +o()K2)| 

If-fll- 


y  y 


00 

/ 


r(x)  F  (x)F  (x)  [  dK  K  J  (KR)  J  (KR) 


exp 


{■[./ 


o  o 


?  9 ~\^/ 2 

4X  +  0  +  (1  +o()K  v/" 


By  using  the  integral  expression  (15),  with  the  source  term  given  by  (10), 

and  integrating  over  the  complete  range  for  x,  ^  and  R  the  following 

o  10- 

expression  is  obtained  for  f  , 


C  N  eEB 
f°  v  o 

4*  rt 


l/2 


4TTD  P  (3/2) 


exp  [(xo  -?>/  ^  ^  e  ° 


Fv  (xo>  F7  <x> 


B 


OKKIo(KR)|-J±0y  +K2|  eXp 


1  -l/Z  f  f  1+4«B  2\ 

J  SXP  *  K  ) 


l/2 


1/2 


In  this  expression  C  ,  which  is  the  normalization  constant  for  the 
maxwell  type  distribution,  exp  C  - (B  6  /(I  +O0)|  ,  is  Cv=  (tob/ 217(1  +<*)) 
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and  is  the  diffusion  coefficient  that  corresponds  to  this  distribution, 
i.  e.  ,  =  (1  +5()^4mB.  The  convenience  in  expressing  the  constant 

in  this  way  will  be  apparent  when  the  far-distant  part  or  lowest  mode 
of  this  complete  distribution  is  discussed  below. 

The  integration  over  K  can  be  carried  out  immediately  by 
using  the  following  equation  given  by  Watson^ 


-l/2 


dKK  Jq(KR)  (K2  +  82)  '  exp  |  -  (K2  +  02) 


■[- 


l/2 


ixi] 


-1/2 


=  (X2  +  R2)  '  exp  -  0  (X2  +  R2)l/2 


The  resulting  expression  for  f  is" 


C  N  eEB 

„0  VO 


1/2 


4TfD  ^(3/2) 


(±+2L)  ri±0L 

\*  1  L°< 


+  u  >A  2  2 

(^-Xo)  +R 


-l/2 


B.  Lowest  Mode 

The  lowest  mode  of  the  distribution  corresponds  to  the 
eigenvalue  (3  =  0(y  =1)  and  to  the  energy  function 


F  (x)  =  e  L 
o  o 


-x  ,  (1/2),  V  /..  -x 


(x)/nq  =  e"X  P {3f2)  j 


N 


This  part  of  the  distribution  is 


When  this  is  expressed  in  terms  of  the  variables  €  ,  z  and  p  it  becomes 


It  will  be  shown  below  that  the  higher  modes  decay  exponentially  with 
distance  from  the  source.  Therefore,  the  lowest  mode  is  the  only  part 
of  the  complete  distribution  that  is  of  interest  at  far  distance  from  the 
source.  To  contrast  the  behavior  of  this  far-distant  part  with  the  usual 
distribution  as  discussed  in  the  previous  section,  it  is  best  to  go  to 


the  limiting  cases  of  low  and  high  field. 
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Low  Field  Limit 


For  low  fields  <K>i  and  f  °  from  (18)  becomes 


>\  -  Cv  -€/kT 

[7  £-eo\\  2I 

J*>r  4Vdv 

^  kTeEB  y  p  J 

-l/2 


exp 


[W,  - 

2kT\  kTeEB/ 


exp* 


eE 

2kT 


rv  ^%\2+  2I1/ 

kTeEB  J  p  J 


If  the  position  is  far  from  the  source,  then 
eEBz  (€  -  £  ) /k 


'kT 


In  this  limit  the  distribution  function  agrees  with  the  usual  distribution 
where  the  density  is  given  by  (8)  with  d/|u  =  kT/e.  Of  course  with 
E  =  0,  the  diffusion  limit  results,  i.  e.  , 


°)  =  (C  ■  Alfo  )(z2  +  p2) 

/(X->oo  v'  v 


1/2 


exp 


(-  6/k T) 


Therefore,  as  was  shown  by  the  qualitative  discussion  above,  when 
the  electrons  are  in  equilibrium  with  the  gas  the  effect  of  gradient 


terms  can  be  neglected. 


i 


1 


I 


! 

i 
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High  Field  Limit 


In  this  case  1  and  f  from  (18)  becomes 

o 


,  o 

i 

o  , 


^  C^eE 


-b6 


Qg  -  €Tq  -  eEz)2  +<X(eEp): 


-1/2 


.  .  exp  -(B/ 2<X)(6-  6  -  eEz) 

o 


[■ 


(€  -  6'o  -  eEz)2  +  0C(eEp)2J  ^  ^ 


(B/20 Q 


Now 


Q€  -  -  eEz)2  +0l(eEp)j  =■  j<£  -  (~Q  -  eEz 


and  therefore  f  can  be  written  as 
o 


/<*<! 


eE  _RC  l 

«5T  e  le_  £o  - eEz 


exp 


^-(B/200  j^(€  -  6q  -  eEz)  +  |.€  -  6q  -  eEz  |  ^ 


exp 


jjB  (eEp)2/4  \€  -  6q  -  eEz  Q  , 


(19) 


(20) 


where  only  terms  that  are  independent  of  <X  or  vary  as  0(.  are  retained. 


It  should  be  apparent  from  (20)  that  the  fraction  of  electrons  with 
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e>e0  +  eEz  becomes  smaller  as G(, decreases  while  the  fraction  with 

£  <  6  +  eEz  is  independent  of  <%.  Therefore  the  limiting  form  for  f  ° 
o  o 

when  Od  is  small  can  be  expressed  as 


f  °)  = 
°  A<! 


Cm  e 


-B6 


4TtDV  Z 


e  -  e 

_ o 

eEz 


(21) 


when  £  <  €  +  eEz,  and 
o 


when  £  >  CQ  +  eEz. 

This  distribution  must  now  be  contrasted  with  usual  distribution 
which,  in  the  same  limit,  is 


(22) 


where  the  density  factor  is  given  by  (8)  with  D/\±  =  1/eB. 

It  is  clear  that  since  (21)  cannot  be  written  as  a  function  of 
energy  times  a  function  of  position,  the  distribution  in  energy  will 
depend  on  position.  To  illustrate  this  new  behavior  the  ratio  of  (21) 
to  (22)  is  plotted  in  Fig.  1  as  function  of  £/eEz  for  various  values  of 
P/  z.  In  this  plot  the  parameter  3/2  eEzB,  which  is  shown  below  to  be 
equal  to  S'.  /eEz  for  on-axis  points  far  from  the  source,  is  taken  to  be 

aV' 

0.1.  Also,  for  simplicity,  the  initial  energy  P  is  taken  equal  to  zero. 

o 
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This  plot  indicates  that  electrons  close  to  the  axis  have  a  higher 
average  energy  than  would  be  expected  while  those  sufficiently  off 
axis  have  a  lower  average  energy.  Also  there  are  fewer  electrons  off 
axis  than  would  have  been  predicted  by  the  usual  distribution. 

If  the  position  is  far  from  the  source,  then  the  quantity 
£/eEz  can  be  considered  small  in  comparispn  to  unity  and  taking  up 
to  first  order  in  this  quantity,  the  distribution  from  (21),  with  =  0, 
becomes 


2 

We  immediately  see  that  when  €/eEz  and  (P/z)  can  be  neglected  in 
comparison  to  unity,  i.  e. ,  for  positions  far  from  the  source  but  very 
close  to  the  axis,  the  distribution  agrees  with  (22).  This  agreement  is 

t 

consistent  with  the  earlier  qualitative  considerations.  The  above  expres¬ 
sion  can  now  be  used  to  obtain  approximate  expressions  for  the  average 

12 

energy  and  density.  These  are 


3(l  +  -l/eEBz) 

2B  (1+  (p/2z)2) 


(23) 
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and 

(1  +  3/2  eEBz) 

n  ~  - - — - - 

4^2(1  +  (p/2z)2) 

These  express,  in  a  quantitative  way,  the  behavior  displayed,  in  Fig.  1. 

Again  we  see  that  the  average  energy  for  electrons  on  axis  is  larger 

then  would  be  predicted  by  the  usual  distribution,  and  by  a  factor  of 

(1  +  1/eEBz).  This  points  out  the  fact,  which  was  not  clear  from  Fig.  1, 

that  for  positions  far  distant  from  the  source  the  average  energy  of  on- 

axis  electrons  goes  to  the  expected  value  of  3/2  B.  Also  we  see  that 

electrons  sufficiently  off  the  axis  have  a  lower  energy  than  expected. 

For  example,  at  an  angle  of  45°  (P/z  =  1)  the  average  energy  is  down 

by  25%,  when  the  term  1/eEBz  is  neglected.  The  usual  expression  for 

the  density  around  a  point  source,  as  contained  in  (22),  can  be  expanded 
2 

in  powers  of  (p/z)  to  be  compared  most  easily  with  the  above  "exact" 
expression  for  the  density.  Such  an  expansion  results  in 


") 


approx. 


-.1-1/2  (P/z)2  +  (1/16)  eEBz  (P/z)4 
4VDv  z 


exp 


eEBP 

4z 


2  1 


and  the  ratio  of  the  densities  can  then  be  expressed  as 


n) 


"  exact" . 
approx. 


4 


~('l  +  1/8  (p/z)2  -  (1/16)  eEBz  (P/z)4  +  3/2  eEBz 


) 
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This  relationship  shows  that  the  density  for  on-axis  electrons  can  be 
higher  than  predicted  by  the  usual  theory,  while  for  electrons  suffi¬ 
ciently  off  axis  the  density  can  be  lower  than  expected.  This  is  in 
agreement  with  the  behavior  shown  in  Fig.  1. 

The  above  results  would  indicate  that  the  average  electron 
energy  would  continue  to  decrease  without  limit  as  the  position  got 
further  and  further  off  the  geometrical  axis.  The  apparent  lack  of  a 
lower  limit  to  the  average  energy  arises  because  the  results  were  obtained 
in  the  limit  of  very  small  By  returning  to  the  expression  for  the 
lowest  mode  as  given  in  (19),  it  can  be  shown  that  for  large  enough p, 
such  that  0(p»  z  +  '(€:  -E  )/eE  i  ,  the  limiting  mean  energy  is  3kT. 

A  qualitative  explanation  as  to  why  the  "exact"  theory 
predicts  an  average  energy  that  can.  differ  from  the  usual  position  inde¬ 
pendent  value  of  3/2B  can  be  given  as  follows.  The  current  density  at 
a  given  point,  which  is- made  up  of  the  diffusion  current  plus  the  drift 
current,  is  a  direct  measure  of  how  asymmetric  the  distribution  is,  i.  e. , 
a  measure  of  how  many  more  electrons  are  moving  in  the  direction  of 
the  current  than  against  it.  Now  the  electrons,  as  a  whole,  gain  energy 
from  the  field  only  because  there  are  more  electrons  moving  against 
the  field  than  with  it,  In  the  usual  theory  in  setting  up  the  balance 
between  the  electrons  gaining  energy  from  the  field  and  losing  energy 
from  collisions  with  the  gas  atoms  it  is  assumed  that  the  contribution 
to  the  electrons  gaining  energy  from  that  part  of  the  asymmetry  corresponding 
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to  the  diffusion  current  can  be  neglected.  Therefore  in  the  "exact" 
theory  where  the  effect  of  the  diffusion  part  of  f^  is  taken  into  account, 
it  is  clear  that  when  the  diffusion  current  adds  to  the  drift  current  the 
mean  energy  will  be  higher  and  when  it  subtracts  from  the  drift  current  • 
the  mean  energy  will  be  lower.  These  conclusions  are  consistent 
with  the  results  for  the  point  source  geometry,  for  in  this  case,  along 
the  axis  the  diffusion  current  aids  the  drift  current  and  it  was  here  that 
the  energy  was  found  to  be  higher,  while  sufficiently  off  the  axis  where 
the  diffusion  current  opposes  the  drift  current  the  energy  was  found  to  be 
lower. 


C.  Higher  Modes 


A  higher  mode  from  (17)  ^in  the  high  field  limit  and  for 
xq  (£  <£q  +  eEz),  is 


f 

V 


C^eEB  N  2 
v _ o 

Air  dvV2(3/2) 


exp 


[V/4  <*o  -?rj 


*  Vo>y*)(*o  -v 

expj^(3(x.-^)J 


When  this  is  written  in  terms  of  the  variables  £  ,  z  and  p  ,  it  becomes 


£  °) 

Voui 


C«eEN  2 

> _ o _ 

4lfDN  p2  (3/2) 


B£ 

e  F  (BS  )F  (B£)(eEz-€  +  6  ) 
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It  is  evident  that  the  higher  modes  have  a  decaying  exponential  in  z. 
From  the  eigenvalue  spectrum  given  above  for  y  ,  it  is  seen  that  the 
characteristic  distance  for  the  first  mode  above  the  fundamental  is 


1 

31  eEB 


=  4.  9/eEB 


Therefore  a  distance  of  the  order  of  1/eEB  must  be  reached  before  the 
lowest  mode  becomes  the  most  important  term.  This  explains  why  the 
discussion  of  the  lowest  mode  as  representing  the  far-distant  distri¬ 
bution  was  carried  out  for  z  >  1/eEB. 


IV.  CONSTANT  CROSS  SECTION 


.  6 

gas  is 


The  equation  that  f° 


must  satisfy  for  a  constant  cross  section 
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+  €  M 
6m 
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v  2  32f° 
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mX2 
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mXsIl  £J 

2  (2m  €  )1//2 


(25) 


with  X=  1/N(T,  where  N  is  the  gas  density  and  C" ,  is  the  momentum 
transfer  cross  section.  This  equation  could  not  be  converted  to  a 
separable  form  by  a  change  to  new  independent  variables  as  was  possible 
in  the  constant  V  case.  However,  when  the  equation  is  taken  to  the 
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high  field  limit,  then  such  a  conversion  can  be  made  and,  of  course, 
in  the  zero  field  limit  it  is  directly  separable.  These  two  limits  will 
therefore  be  presented  as  separate  problems. 


1  c) 

6  3€ 


A.  Zero  Field  Limit 


In  this  limit  the  equation  for  f  is 


[e 2  (f° 


dfOV  ! 

+  ^  96 


'J 


+  V2£° 

6m  OL 


Xms 

2(2mS  )1//2 


This  is  to  be  solved  with  the  source  function,  S,  given  by  (10).  The 
complete  distribution  function  for  this  case  is  obtained  by  the  same 
procedure  as  was  used  with  constant  V  and  is 


f°  = 


41,D°o  5=0 


00  €r/kT  2  2 
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(26) 


The  energy  eigenfunction,  Fj^  (u),  are  solutions  to  the  separated  energy 
equation 

% 

U^T  +  (2  +  u)  (2  +JL  )  F  =  0  , 

,  2  du 

du 


where  u  =  6  /kT.  These  functions  are  discrete  and  given  by 
F^  (u)  =  e"U  (1)  (u)/N£ 
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j  * 


i 


with  Si  =  0,  1,  2  .  .  .  .  ao.  N^,.  which  is  the  normalization  constant 
for  the  density  function  r(u)  =  ueU,  is  given  by 

Njj, 2  .=  (i  +4)  p2  U+-2). 


The  constant  in  front  is  expressed  in  terms  of  C  ,  the  normalization 

o 


constant  for  the  maxwell  distribution  e 


C/kT 


and  in  terms  of  Dq,  the  corresponding  diffusion  coefficient,  given 
by  Dq  -(2  kT/m  if)^2.  The  lowest  mode  from  (26)  is  of  the 

expected  form 


(Co/4^Dq)  e  e/kT(p2  +  z2) 

The  higher  modes  decay  exponentially  with  distance  from  the  source 
with  the  characteristic  length  of 


B.  High  Field  Limit 

In  Eq.  (25),  for  the  limit  of  high  field,  the  term  representing 
electrons  gaining  energy  from  collisions  with  the  gas  atoms  can  be 
neglected  and  the  resulting  equation  is 
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It  is  clear  that  because  of  the  mixed  derivatives  in  £  and  z  the  equation 
cannot  be  separated.  However,  by  changing  to  new  dimensionless 
variables  given  by 

y  =  (A£  ) 2  , 

7  =  2A(  £  -  eEz)  , 
and  7  =  2  eEA  p  , 

1/2 

where  A  =  (1/eEX  )(3m/M)  '  ,  the  equation  becomes 
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This  is  now  in  a  separable  form.  The  formal  solution  to  this  equation, 
which  can  be  obtained  by  the  same  procedure  that  was  used  for  constant  - 
collision  frequency,  is 


KJo(Klj) 


eS 


i  <Kk,Fj>  (r2-  y> 


1/2  •  1/2 
for  2  yQ  and  equal  to  zero  for  2  yQ  ^  .  The  constant  Cg-  is 

the  usual  normalization  constant  for  the  Druyvesteyn  energy  function, 

exp  [-(as)2:  ,  and  is 
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and  D  is  the  corresponding  diffusion  coefficient 
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The  functions  Fj^  (K  ,  y)  are  the  eigenfunctions  of  the  separated  energy 
equation. 


d2Fo  dF«  9  \/9 

+  (i  +  y)T£-  +  d-K2  +  yi/^)FJ,  =0  , 


dy 


dy 


«  2 

with  corresponding  eigenvalues  (K ).  The  density  function  for  this 


+  .  .  1/2  y 

equation  is  y  e  . 


Because  the  behavior  of  this  equation  at  its  irregular  singular 

13 

point  at  y  =  oo  is  of  a  more  complicated  type  than  the.  usual  exponential, 
the  exact  analytic  form  of  those  solution  could  not  be  obtained.  However, 
an  approximate  expression  for  the  lowest  eigenfunction  can  be  obtained 
by  a  perturbation  calculation,  in  a  form  which  is  appropriate  for  expressing 


the  far-distant  behavior  of  the  distribution  function. 
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The  far-distant  mode  of  the  complete  distribution  is  given  by 


CT  eEA 

'TT  Dr 


°o 

°  pKKJ0(KH)E'6(fdy0)Fa(£y): 


We  will  consider  that  Fq  and  can  be  expanded  into  a  power  series 

.  „2  , 

in  K  ,  i.  e. , 


F  (K2,  y)  =  F  *°)  +  K2F  ^  +  K4F  ^  +  . 

O  O’’  o  o 

and 

(K2)  =  +  K2^^  +  K4  + 


When  these  are  inserted  into  the  above  integral,  fQ  becomes 


f  °- 


C^-  eEA  y 


o  TTD, 


oo 

°  /  dK  K  Jo  (KH(Fo(°)(>'o)+K2Fo(1>(yo)+-  •  ) 


(F  <°>(y)  ♦  K2  F  «M  ♦  .  .  .  .)  expr-K2(2yo1/2  -? )  If' "1 


x  exp 


(2yo1/2-  ?)  < 1  -  K4(2yo1/2  -  <?)  t Po(2)  + 

(l-K6(2yo1/2-^)^fo(3)H- 


«  •  0 


These  series  can  now  be  multiplied  out  and  f  can  then  be  expressed  as 
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Cn.eE A  y  ~  .*77  9* 

£o' ’  e  0  <*KKJo(Kf(g( 


(1  +  K2g1  +  K4g2  +  .  .  .  )  exp  ( 2 y^1'  2  -^)  K2  tpj^ 


(27) 


where  the  g-functions  are  given  by 


g  =  F  ^o)(y  )  F  ^(y) 
yo  o  o  o  * 


9‘  *  Fo(\> +  Fo<o)(y> 


Fo(2)(y)  fo(2)(y0) 


Fo(l)(y)  Fo(l)(yo) 

F  '~'(y  )  F  (o)(y)  F  (o)(y  ) 
o  o  o  o  o 


^/S-i)€\etc. 


Eq.  (27)  Can  then  be  integrated  term  by  term  using  the  following  bessel. 
integral  formula  given  by  Watson,  ^ 


The  expression  for  f  °  becomes 
o 
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_  O  C<T  sea/0  exp  [-  (2yo1/2  -'J )  exp  \j-y\M2yV2  -<•) 

£°  "  2Wlv  (2y j/2.^)  I f)V) 

(  giLi(p)(n2/4(2yn1/2-?)lPoaj)  92L2(0)(~n.2/il2yo1^  -%)  VJ0)  ~) 

t  WWf  +"'y 


If  the  region  of  interest  is  restricted  to  points  that  are  far  distant  from 
the  source  but  close  to  the  axis,  then  the  quantities 


(p/z)2,  6/eEz,  and  1/eEAz 


can  be  considered  small  in  comparison  to  unity.  When  the  Laguerre 

1/2  -M9 

polynomials  are  expanded  in  (28)  and  it  is  recalled  that  (2yQ  --sj) 
when  expressed  in  terms  of  £  and  z  is  given  by 


(2yo^2  -  ^)  =  2  eEzA^l  -  (€ 


then  the  terms  that  are  first  order  in  the  above  small  quantities  can  be 
picked  out.  Therefore,  the  most  important  terms  in  the  bracket  are 
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Thus  knowing  Fq  to  first  order  in  K  and  to  second  order  in  K 

results  in  an  expression  for  f°  that  is  valid  near  the  axis  (first  order 
2 

in  (P/z)  }and  at  far  distances  from  the  source  (first  order  in  1/eEAz). 

In  the  appendix,  Fq  and  (^q  are  obtained  by  a  perturbation 

2 

calculation  up  to  first  and  second  order  in  K  ,  respectively.  From 
these  results  the  g-functions  and  the  expansion  coefficients  for  ip 
can  be  obtained  as 


and 


» 


Values  for  a^,  for  n  up  to  four,  are  given  in  Table  II. of  the  appendix. 

Therefore  fQ0  can  be  expressed  in  terms  of  &  ,  z  and  p  as 


CL-exp  -(A6) 

f  =  - ;(l-(C-€  VeEz)  exp 


^TDg.  z 


ir1//2  eEA/^ 

4  z  (j-(€  -  €Q)/eEz  )  j 


3/2  1/2 

av  vjreEAz  .  J  M  g, 


■id  Tf  g,  \  a  ir  etAz  .  ,  /  a 


eEzA  \  2  1/2)  *  * 

(29) 


x 
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It  is  immediately  obvious  that  for  points  very  near  the  axis  but  at  far 
distance  from  the  source,  the  expression  reduces  to  the  usual  approxi¬ 
mate  form  taken  to  this  same  limit,  i.  e. , 

(%  ATTD^z)  exp  ^-(Ae  )^  exp  J 

1/2 

where  ~  ‘]y  eA. 

For  points  off  the  axis  it  is  more  difficult  in  this  case  to 
compare  (29)  with  the  approximate  f°  than  it  was  for  the  case  of  constant 
collision  frequency.  However,  expressions  can  be  obtained  for  the 
average  energy,  £  ,  and  for  the  density,  n,  which  can  be  compared 

3V 

with  those  from  the  approximate  f°.  These  expressions  are 


£  -  P(5/4)  j 

av  ~  P(  3/4)  A 


^-0.198  (P/*)2-^"] 


and 


n  *  (WD^z) 


x  exp 


-1 


1  -  0.204  (P/z)2  +  "^f|+  .02.3  eEAz  (P/z)4 


lVl//2  eEA P2 
4z 


i  ” 


The  integrals  involving  the  function  were  evaluated  by  using  the 

g 

generating  function  for  the  Laguerre  polynomials.  The  numerical 
constants  in  these  expressions  are  accurate  to  a  few  per  cent.  As  in 
the  constant  collision  frequency  case  the  initial  energy,  0  ,  was 
taken  to  be  zero.  It  is  to  be  noted  that  £  goes  to  the  expected 
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value  of  (P( 5/4)/P(3/4)  Ay  or  points  very  close  to  the  axis  but  far 
from  the  source.  Also  it  is  apparent  from  the  comparison  of  these 
relationships  with  those  obtained  for  constant  collision  frequency, 
which  are  given  by  (23)  and  (24),  that  the  qualitative  discussion  and 
conclusions  for  the  case  of  constant  collision  frequency  will  apply 
directly  to  the  constant  cross  section  case. 


VI.  TOWNSEND  TYPE  D/(i  EXPERIMENT 

The  theory  as  developed  above  has  demonstrated  that  in  general 
it  is  not  correct  to  ascribe  to  electrons  moving  through  a  gas  a  distri¬ 
bution  in  energy  that  is  independent  of  position.  In  turn  this  implies 
that  such  electrons  are  not  characterized  by  a  unique  diffusion  coeffi¬ 
cient  and  mobility.  However,  the  assumption  of  a  unique  D  and  ^  does 

form  the  basis  for  the  usual  interpretation  of  the  Townsend  type  D/[i 
1'  2 

experiment.  1  Therefore,  this  type  of  experiment  will  be  re-examined 
in  terms  of  the  above  theory  to  find  if,  under  the  conditions  of  the  actual 
experiments,  appreciable  errors  are  introduced  by  using  the  usual  inter¬ 
pretation.  While  it  is  recognized  that  the  above  theory  does  not  take 
into  account  the  effect  of  electrode  boundaries  or  of  inelastic  processes 
in  the  gas,  both  of  which  are  important  for  an  accurate  description  of  the 
experiments  as  actually  carried  out,  still  the  essential  features  of  the 
D/fj.  experiment  are  represented  by  the  point  source  in  a  gas  in  which 
only  elastic  collisions  take  place.  It  is  therefore  reasonable  to  suppose 
that  the  theory  will  give  at  least  an  estimate  of  the  errors  involved. 
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When  the  distribution  in  energy  is  assumed  to  be  independent 
of  position  for  a  constant  collision  frequency  gas,  the  expression  for 
the  density  normalized  with  respect  to  the  value  on  axis  can  be  obtained 
from  (8).  This  density  ratio,  which  will  be  denoted  by  N,  is  given  by 

^Approx.  ■{*  +  <  P/z>2) 

This  ratio  could  just  as  well  be  expressed  in  terms  of  D/p,  where 

D/p.  =  1/eB,  This  expression  is  taken  for  the  high  field  limit  since 

it  has  been  shown  earlier  that  the  usual  expression  for  the  density  is 

only  in  error  when  the  average  electron  energy  is  large  in  comparison 

to  thermal  energies.  It  is  clear  that  if  the  ratio  of  the  density  at  some 

point  off  axis  to  the  density  on  axis  is  measured,  i.  e. ,  if  the  density 

14 

ratio  is  measured,  then  B  (or  D/u)  can  be  calculated. 

The  more  accurate  description  of  this  experimental  situation, 
has  shown  that  while  the  density  ratio  is  still  given  by  the  above  expres¬ 
sion  for  points  very  close  to  the  axis,  a  more  correct  expression  for 
points  further  off  the  axis  is 

N) exact  =  [/  1/4  (P/Z>2]  exp  [j  ^  (P/z)D  • 

These  expressions  can  most  easily  be  contrasted  by  expanding  each  in 
2 

powers  of  (P/z)  up  to  the  point  that  their  difference  is  apparent.  This 


(l+  (P/z)2) 


■1/2 


results  in 


39 


N)aPProx.  '  (j  -  ‘/2  (PA)4)exp  ^  (P/z)2] 


(30) 


and  in 


^ exact  =  (*  ■  3/8  (P/z)2)  exp  [j  (P/z)23  • 

(  ¥)■ 


(31) 


Typical  experimental  values  for  (p/z)  and  eEBzl 


which  are 


also  representative  of  previous  measurements,  can  be  obtained  from 

2 

the  recent  work  of  R.  W.  Warren  and  the  author.  In  these  measurements 
the  two  experimental  tubes  that  were  used  differed  in  z,  the  distance 
from  the  source  to  the  plane  of  measurement,  and  in  (P/z).  Table  I 
gives  the  range  of  p.Ez/D  covered  in  these  measurements  along  with 
the  values  of  (P/z)  (two  values  per  tube)  for  each  of  the  two  tubes. 

Also  given  in  Table  I  is  the  maximum  fractional  difference  in  N,  i.  e. , 


rM  n)  -  N ) 

6  N  _  '  agorae.  S 

' 


exact 


N 


ipprox. 


2 

as  calculated  from  (30)  and  (31)  for  each  (^/z)  .  These  figures  show 
that  for  a  given  value  of  eEBz  the  density  ratio  as  predicted  by  the  two 
different  theories  can  differ  appreciably  for  typical  experimental  condi¬ 
tions. 

However,  the  important  question  is  how  different  is  the  value 
of  p. E/D  as  predicted  by  the  two  theories  for  a  given  N.  An  expression 


:40 


for  the  fractional  difference  in  D/\x  (or  ~r)  for  a  given  N  can  be  obtained 

eb 

by  using  (30)  and  (31).  This  is  given  by 


This  can  be  expressed  as 


S(D/u)  „  1  D 

D/p  2  pEz 


2 

where  In  (1/N)  has  been  replaced  by  (pEz/4D)(  p/z)  .  Table  I  gives  the 

maximum  value  that  this  fractional  difference  can  attain.  It  is  clear 

from  these  figures  that  while  the  fractional  difference  in  D/p  is  not  as 

large  as  in  N,  still  the  difference  of  the  order  20%  corresponding  to 

2 

the  largest  value  of  <P /*>  should  be  experimentally  observable.  In 
the  course  of  the  measurements  described  in  Ref.  2,  certain  inconsis¬ 
tencies  did  arise  using  the  conventional  interpretation  of  this  experiment. 
However.-  these  inconsistencies,  which  were  eventually  resolved  by  an 
empirical  approach,  could  not  be  explained,  even  qualitatively,  by  the 
results  of  the  present  theoretical  investigation.  Therefore  it  would 
appear  that  in  most  cases  the  experiments  have  not  been  appreciably 

affected  by  using  the  usual  interpretation  of  this  experiment  and  in  the 
2 

cases  where  ( p/z  r  was  large  enough  for  appreciable  deviations  to 
exist,  such  deviations  were  masked  by  other  effects. 
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APPENDIX 

Lowest  Mode  Energy  Function  for  Constant  'T“ 

The  energy  differential  equation  for  this  case,  is 

y^-T+  (1  +  y)  +  (1  -  K2  +  y1/2(p )  F  =  0. 
dy^ '  ay 

Making  the  transformation 

F(y)  =  e"V  h(y)  ’ 

the  equation  becomes 

Y  +  (1  "  y)  djj  +  (yl/2  ^  “  r2)  h  =  0  •  (A- 1} 

We  want  to  obtain  the  first  few  terms  of  a  power  series  expansion  in 
2 

K  for  tfo  and  hQ,  the  lowest  eigenvalue  and  eigenfunction  for  this 
equation.  That  is,  we  want  to  find  the  first  few  terms  in 

h 

o 

and  in  ^ 


-  h  (0)  +  K2  h  (1)  +  K4  h  (2>  +  . 
o  o  o 


•  ) 


=  <fj°>  +  K2^')  +  K4  lfo<2)  + . 


A-l 


A3 


—  Here  we  have  considered  to  be  the  eigenvalue  which 
2 

is  a  function  of  K  .  However,  in  the  actual  calculation  of  above 

i 

\ 

terms  it  is  more  convenient  to  reverse  this  viewpoint  and  to  consider 
2 

K  the  eigenvalue  which  is  a  function  of  ^  and  then  to  consider 
1/2 

the  term(y.  )  in  (A„  1)  as  a  perturbation.  It  is  convenient  to  put 
(A.  1)  into  a  standard  quantum  mechanical  form^5  so  that  the  usual 
perturbation  formulas  can  be  used.  Eq.  (A.  1)  can  then  be  expressed 
as 


H  11J  =w  Vi) 

o  To  '  To  oTo 


where 


dy 


H1  =  y 


1/2 


and  w  =  K 
o 

When  IP  is  equal  to  zero,  the  complete  set  of  discrete,  orthonormal 
eigenfunctions  are^ 

Un  =  Ln(o):(y)/nl  , 

where  L  ^  is  a  Laguerre  polynomial  with  corresponding  eigenvalues  of 
n 


A- 2 
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E  =  0,  -1,  -2.,  ...  -xi  ...  -oo  , 


These  functions  are  orthogonal  with  respect  to  the  density  function 
e'y. 

Using  the  usual  perturbation  formulas,  4^  and  w  ,  to  second 

*  o  o 

15 


order  in  perturbation  theory,  are' 

oo 

w0  =  'Ho "  y22" 


03  (H'  )2 

no 


n^l  n 


oo 


and  U 


H’ 


no 


,  E 
n=l  n 


U 


where  we  have  explicitly  put  Eq  =  o  and  with  H1^  given  by 


oo 

/ 


H’  =  /  y1//2  e"Y  U  U  dy  . 
no  J  on 


The  value  of  H1  can  be  evaluated  with  the  help  of  the  generating 

9 

function  for  the  Laguerre  polynomials  and  this  quantity  is  tabulated 
in  Table-  II  for  n  up  to  four. 

We  can  now  turn  around  and  obtain  (^(w)  and  ^Q(w)  and 
these  relations  are 


Y  = 


W  _ w: 


00  (H*  -  )2 

no 


o  H'ft.  (HJjr)  v 


+  .  . 


"00  Vil0Q7  n=l  n 


A- 3 
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and 


00  H  1 

-na_  „  + 

o  o  Hj.  £—•  £  n 

OQ  n=l  n 


4>  -u  J~ 

*  o  o  H,.  fc— ■ 


The  function  which  is  normalized  with  respect  to  e  y  is  now 

renormalized  to  y^  2  e  y.  Then  U)  and  h  can  be  written  in  terms 

To  o 

of  K2  as 


\  00 

ho  -  K2  %  L^(y) 


and 


*f0  =  (2 /r^2)  K2  (1  -  K2  aQ)  . 


The  constants  a  can  be  expressed  in  terms  of  H'  as 
n  no 


=-  z 

or*  <-  »"■ 


°°  (H'  )2 

no 


o  'IT  n 

n=l 


and 


2  .H'no.  f  .  , 

■ —  for  n  ^  1  . 


3n  ^yl/2  nnl 


Values  for  a^  are  given  in  Table  II  for  n  up  to  four.  In  summary,  the 
expansion  coefficients,  for  hQ  and  are 


A- 4 


TABLE  I.  Maximum  Fractional  Difference  in  n  and  D/|j. 


uEz 

D 


(|°/z)2 


50  -  500  (Long  Tube) 


0.  005 

0.  3 

0.  1 

0.  02 

5 

0.  5 

0.  24 

15 

8 

0.9 

100 

o 

22 

2-10  (Short  Tube) 


TABLE  II.  Constants  for  the  Perturbation  Calculation 


n 

2H'  /fr1''2' 

no  11 

a 

n 

0 

1 

fcS  0. 259 

1 

-1/2 

1/2 

2 

-1/8 

1/32 

3 

-3/48 

1/288 

4 

-15/384 

5/12,  288 
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"Spatially  Dependent  Energy  Distributions  for  Electrons  Drifting  Through  a  Gas 
in  a  Uniform  Electric  field"  by  James  H.  Parker,  Jr.  (Scientific  Paper 
63-928-113-P2) 


Page  12 


In  eq.  (9)  the  term  on  the  left  hand  side  should  read  (-  . 


Page  17 


read 


The  denominator  in  the  first  equation  at  the  top  of  the  page  should 

"  Q  +  4oc  <p  +  (l  +oo  i <f]  1 ^  ". 


Page  18 


sign. 


In  eq.  17  the  function  F  (x  )  should  be  inserted  after  the  summation 

7  0 


Page  20 


exp 


E 


In  line  8  from  the  top  the  exponential  term  should  read 

1/2 


eE  ,_2  .  „2 
2kT 


C*  +  p) 


Page  25 


The  inequality  on  line  11  should  read  "  <xV^»  z  +  -  G^eEj  ", 


Page  26 


The  inequality  in  bracket  on  line  13  should  read  "(or  6  <  £  +  eEz)" 

o 
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Paee  31 

The  first  sentence  of  the  last  paragraph  should  read,  ’’Because  the 
behavior  of  this  equation  near  its  irregular  singular  point  at  y  =  do  requires 
a  more  complicated  solution  than  the  usual  exponential,  the  exact  analytic  form 
of  the  solutions  could  not  be  obtained". 
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r  2  / 

At  the  end  of  reference  12  replace  "1/  j^l  -  (p/2z)  "  by 


[l  -  (P/2z)2J  (l  -  2eEBzV'* 


tt 


